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Abstract. We investigate the behaviour of exact closed bouncing Friedmann 
universes in theories with varying constants. We show that the simplest BSBM 
varying-alpha theory leads to a bouncing universe. The value of alpha increases 
monotonically, remaining approximately constant during most of each cycle, but 
s^ , increasing significantly around each bounce. When dissipation is introduced we show 

Qs ' that in each new cycle the universe expands for longer and to a larger size. We find a 

VO ■ similar effect for closed bouncing universes in Brans-Dicke theory, where G also varies 

~^ I monotonically in time from cycle to cycle. Similar behaviour occurs also in varying 

^—s ■ speed of light theories. 
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The "bouncing universe" is a modern cosmological reincarnation of an ancient 
fascination with the cychc patterns of nature and the myth of the "eternal return" 
PI El 0] • Gravitation theories hke general relativity allow us to make precise models 
of this popular conception of a "phoenix" cosmology, in which a universe periodically 
collapses to a Big Crunch, only to rebound into a new state of expansion, as if emerging 
from a unique Big Bang |1]. Multiple bounces are possible but each cycle lasts longer 
and expands to a larger maximum size than the previous one, a consequence of a simple 
application of the second law of thermodynamics IH El El EI, unless there is a finite 
positive cosmological constant, in which case the oscillations must eventually cease jH] 
and are replaced by eternal de Sitter expansion. A sequence of many oscillations will 
drive the bouncing closed universe closer and closer to flatness. 

"Quantum gravity" effects are invariably invoked to justify the bounce; possible 
detailed calculations, however, have only recently emerged. In loop quantum gravity, 
the semi-classical Friedmann equations receive corrections that produce a bounce PI ITU]. 
The ekpyrotic model of the universe, inspired by string/M-theory, is another possible 
realization of phoenix cosmology JT]. It is also possible that ghost fields - fields 
endowed with negative energy - are capable of producing a classical bounce (this idea 
has been often rediscovered; see [12] for a good review). Classical bounces produced by 
conventional scalar fields with potentials which only violate the strong energy condition 
are difficult to produce in universes that grow large enough to be realistic: typically the 
probability of bounce is of order the ratio of the minimum to the maximum expansion 

size jinnnj. 

It has been speculated that whatever causes a collapsing universe to bounce can 
reprocesses some aspects of physics, either randomly [1^], or systematically [iH], by 
changing the particle spectrum or resetting the dimensionless "constants" of Nature. 
Both of these options are severely constrained by anthropic requirements but it is 
interesting to ask whether there are monotonic or asymptotic trends in the values of 
some quantities, as seems to be the case for the degree of flatness of the universe, over 
many bounces. This matter is clearly of great importance in the context of varying- 
constant theories [T71 ITH] . Here, quantities which are traditionally constants become 
space-time variables and if singularities are avoided in the bounce then their evolution 
from cycle to cycle is predictable by the field equations rather than the outcome of 
effectively random reprocessing. The values of any dimensionless 'constants' of Nature 
could evolve towards asymptotic attractors if they are allowed to be variables in a self- 
consistent theory. Studies of these theories are also important in assessing the stability 
and level of fiuctuations in a bouncing universe. It has been suggested that thermal 
fiuctuations in bouncing models could be the origin of the cosmic structure EH] ; this 
would provide a distinct alternative to an origin from vacuum quantum fiuctuations in 
a de Sitter phase of cosmological expansion. 

Although these claims are intriguing, it is difficult to evaluate them in the absence 
of a concrete model for the bounce, which is usually viewed as a black box from 
which anything can emerge ^3]. In this paper we examine some exactly soluble 



Bouncing Universes with Varying Constants 3 

examples. First, we consider a simple bouncing cosmology following from the simplest 
BSBM varying-constant theory we developed (2^ from Bekenstein's varying-alpha 
model |2I]. We derive the result that for suitable couplings (indeed those favoured 
by observations [22123; see however [21j) the theory leads to a bouncing universe. We 
derive similar behaviour in the Brans-Dicke theory of a varying gravitational constant 
G j23] , and in a class of varying speed of light theories [211 I2Z1 I2H1 1211 (HH] ■ In each 
case we are able to find exact solutions in the absence of dissipation and compute 
the evolutionary trends in the fine structure constant and G from cycle to cycle when 
dissipation occurs in accord with the Second Law. We find the interesting result that 
the varying constants in these theories change monotonically from cycle to cycle when 
the scale factor oscillates: the scalar fields determining the constants in each cycle do 
not oscillate. 

Crucial to our models is the idea that cosmological fields may have a negative 
energy. Such fields are called ghosts J2] , and are far from new, having found widespread 
application in the study of steady-state cosmology [HI], phantom dark matter [S21ISS1, 
and ft-essence j^l]. Bekenstein's theory |2I] is rooted in the use of a real scalar dielectric 
field ip, representing the allowed variation of the electron charge according to e = eoe^, 
where Cq is the present-day value of the electron charge, and so the fine structure 
'constant' evolves with respect to its present-day value, a^, as a = a^e^'^. By redefining 
the electromagnetic gauge field [201, ^/i, as a^ = e^A^, and the electromagnetic field 
tensor, F^,^, as f^^, = e'^F^^, the action may be written as 

S = J d^x^^ [Cg + C^at + C^ + C,^e-^^) , (1) 

where C^ = —"^d^ipd^ip, Cem = —\ffiuf^'^, and Cmat does not depend on ip. The 
gravitational Lagrangian is the usual Cg = j^^R, with R being the curvature scalar. 
If the scalar coupling satisfies u < 0, then ip has a negative kinetic energy term, and is 
thus a ghost field. 

This theory can fit the observational constraints on varying alpha reported 
by 1^ ES] as well as others [201 IHEI- However for this to be possible in the ghost- 
free case with u; > one has to choose a very special type of dark matter, in which 
the magnetostatic energy B^ dominates over the electrostatic energy E"^. Even though 
a dark-matter candidate was found satisfying this condition (superconducting cosmic 
strings), most types of matter, including baryonic matter, are E"^ dominated, for which 
observational data implies tu < 0, so that ^ is a ghost field. 

Ghosts have been criticized on a variety of grounds. Classically they are a source of 
instabilities if coupled to other forms of matter, since they will try to off-load an infinite 
amount of positive energy into them. This is not necessarily cataclysmic if the rate of 
these processes is sufficiently slow. For instance, in steady-state cosmology to negative 
probabilities. At the quantum level, ghosts lead to negative norm states and so negative 
probabilities. The quantum instabilities are also much more severe and are present even 
without direct coupling to matter, for example in runaway particle production via the 
graviton vertex. Hence, at the quantum level ghosts are pathological. 
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However, we know that quantisation of the field ip is pathological even for c<j > 0. 
For one thing, the theory is non-renormalisable. The attitude to ip should therefore be 
similar to that with regards to gravity: "don't quantise" . General relativity is also, at 
face value, non-renormalisable: for instance, the quantum corrections to the relativistic 
precession of the perihelion of Mercury are infinite. This doesn't stop the classical theory 
from being very successful. It may be that the quantisation of ghosts is simply more 
subtle; ghosts have been found as type //* string theories [S^ . 

Non-relativistic matter and the cosmological constant may be neglected near a 
bounce, so let us consider a Friedmann-Robertson- Walker (FRW) universe filled with 
radiation and a dielectric field, ^p. The cosmological equations for BSBM varying-a 
theory ((T)) are 

H' = l{pre-''' + P,)-^,, (2) 

- = -I [^Pr^-^" + 4p^) , (3) 

where we have set SvrG = 1, i7 = aja is the Hubble expansion rate, K is the 3-curvature 
constant, and p^ = ujiIj'^/2. For the scalar field, in the absence of non-relativistic matter, 
we have 

^ + 3Hip = 0. (4) 

so ^Z) oc a~^. We should not ignore the possibility that at high curvatures quantum 
processes may allow the conversion of ip energy into radiation. We take this into account 
by introducing variable p,. = p^e'"^^ and rewriting the conservation equations as: 

p^ + 6Hp^= - s{pr,ij,a), (5) 

Pr + 4:Hpr = s{pr,ij,a). (6) 

In this case, the equation of motion for ip will contain an additional s term which models 
energy transfer between the ip field and the radiation sea in accord with the second law 
of thermodynamics. We shall consider the implications of such a process, in all its 
generality, below, but let us look first at the equations neglecting this coupling function. 

Consider first a model s = bouncing universe which is exactly soluble. Taking 
K = +1 and uj < 0, Eqn. (0) is 

where S and F are positive constants. In terms of conformal time drj = a^^dt, this can 
be integrated to give 



a\7]) = - [F + VF2 - ASsm{2{r] + r]o)}\ (8) 

when F^ > 4S. Identifying the expansion maximum and minimum, we see that a{ri) is 
given by 



2 1 

" =2 



amax + «min + («max-«min)sin{2(r7 + r7o)} , (9) 
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Figure 1. Examples of solutions for s = 0: a stable, static universe with a increasing 
exponentially (solid line); a stable perturbation around the static solution (dotted), 
and a bouncing universe (dashed). 



where amax is global expansion maximum and a^^^ is the global minimum of a{ri), defined 
by 



2 

max 
2 
min 



r ± vr2 - 4S 



(10) 



Since ip = Ca ^ we have, for u < 0, that S = —ujC"^ /2 and the scalar field driving 
time-variation of the fine structure constant is given by 



i, 



±- — -tan 

\lo\ 



rtan(r/ + r/o) + VP-4^' 

27^ 



In Fig. [T] we plot these solutions and show a{t) and a{t) (recall that a oc e^'^), as 
functions of proper time, t. We note the steady increase of a with time despite the 
oscillatory behaviour of the expansion scale factor. 

When s = we have a variety of oscillating solutions whose characteristics depend 
on the initial conditions. For bouncing solutions a remains nearly constant during 
each cycle but changes sharply, but still monotonically, at the bounce. There is no 
significant change of behaviour at the expansion maximum which also implies that there 
should be no gross difference in evolution inside and outside spherical overdensities far 
from the bounce. With a^ax ^ O'mim and setting F = p^a^^ /?> and S = —p^a^/3, 

= ^/T. We can then see that the bounce duration 



we have a. 
is At 



mm 
2 



2 I 



o/i and ttmax 

^. Since ^ ~ y6F3/2/(5 



u 



l/2x 



near the bounce, we find that 



A^ ~ i/6/ \uj\, independently of initial conditions, during each bounce. 
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Figure 2. A bouncing universe with uj = —0.01 and s = lOQp^. Otlier couplings lead 
to qualitatively identical scenarios. The growth in a{t) and total period from cycle to 
cycle is accompanied by increase in a which occurs primarily in the vicinity of each 
bounce. 



The extreme case is a stable static universe. Setting a = and a = 0, we can see 

-Pr/2, giving a -- 



that this case is reahzed when p^ = —pr/2, giving a = jQ/pr- For such a universe ^p 
evolves linearly in t, and since we have a oc e"^"^ there is exponentially rapid increase 
P7j . Even though such a universe is static, the rulers and clocks of observers change 
as alpha changes, so that they actually observe a Milne universe. We can see that the 
solution is stable because homogeneous and isotropic perturbations lead to a universe 
with regular sinusoidal oscillations as can be seen in Fig. ^ Such solutions are described 
by © in the case where Qmax ~ cimin- This situation differs from that found in general 
relativity in the absence of ghost fields |HHj . 

If s is a non-vanishing then, regardless of its exact functional form, there are two 
type of solutions. If s 7^ at all times, then sooner or later the universe enters a steady- 
state evolution with exponential expansion and constant overall energy density ensured 
by the appropriate transfer of energy between the ip field and radiation. However, we 
expect that these energy-transport processes will switch off at low curvatures, when 
the universe expands to a sufficiently large size (a >> amin) and transport processes 
become collisionless and far slower than the expansion rate. Then, the typical evolution 
is as plotted in Fig. |2l Again, ip is approximately constant during each cycle and 
changes dramatically at the bounce. In addition, each cycle is now bigger than the 
previous one, because F increases at each bounce. This is an interesting realisation of 
the standard Tolman scenario. Cycles get bigger (and entropy is generated near the 
bounce) specifically because radiation is produced from the scalar field close to each 
bounce. In producing Fig. |21we have used s oc p^ {ip decays into radiation), but other 
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Figure 3. A bouncing universe for a Brans-Dicke theory with ujbd < —3/2, in the 
Jordan frame. In this case the field (j) and matter are decoupled: they never exchange 
energy. The amplitude of oscillations in the scale factor grows monotonically in time 
and G increases in the Jordan frame. The increases in G occur primarily near the 
bounce. 



functional forms may be used with similar effects. 

We now examine similar solutions for the Brans-Dicke (BD) theory of varying 
G |2S]- In the Einstein frame, if the matter content is pure radiation, we recover 
the same equations. All we need to do, then, is convert the above results into the 
Jordan frame: the results found for a, p, and ip should then be translated into variables 
a J = a/a/^, Pj = p0^, and (p = e^ (the latter corresponding roughly to l/G). Under 
this transformation we obtain u = ojbd + fj where ojbd in the Brand-Dicke coupling 
parameter. Thus we need ujbd < —3/2 for the Brans-Dicke field to behave like a proper 
ghost in the Einstein frame. The resulting dynamics is plotted in Fig. El 

These results may be understood analytically. The essential BD field equations in 
the Jordan frame are 



H 



+ 3H(j) = 

2 _ 8^P 



-H^ + 



wbd<, 



K 



where overdots now refer to derivatives with respect to tj and H = hj/aj. Hence 
A 



:ii) 



A 



a J Ojcj) a J (j) 



aj0 Wbd<, 



K 



with A and A > constants. We are interested in negative uj solutions with K = 1, 
which give oscillating closed universes (note that wbd < is not sufficient for an 
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expansion minimum, we need ubd < ~3/2). Following the techniques of ^| we put 

C > 0, we have simple exact solution 



^(2cc; + 3) = —C, and in the bouncing case A^ 
in terms of conformal Jordan time: 

c^a] = ^ + l^^ 



C sm{2{r] + r]o)}. (12) 

We see the same behaviour as displayed by the radiation-scalar universe given above 
with s = 0. The minimum value of 0aj is A — VA^ — C and the maximum is A+VA^ — C. 
So in conformal time (pa'j ^ a'j/G undergoes oscillations of increasing amplitude as the 
entropy increases (ie if A increases in value to model increasing radiation density); that 
is, the horizon area ('entropy') in Planck units increases from cycle to cycle. The full 
solution for a{ri) and (j){ri) is then obtained using eqns. ()11|) and p2|) : 

2A 



(p = (pi exp 
with 01 constant and 



C 



arctan 



Atan{?7 + r/o} + VA^ - C 

7c , 



(13) 



a 



J 



^ + lyp3^sin{2(ry + r/o)} 



X 



exp 



2A 



C 



arctan 



A tan{?7 + r/o} + VA^ - C 



C 



(14) 



The effect of increasing radiation entropy can be seen by increasing the constant A in 
these expressions. In this model the field (p and radiation are fully decoupled. Although 
the size of the universe increases in each successive cycle, its size with respect to Planck 
units remains the same, unless of course we consider a model in which the field (p and 
radiation may exchange energy. 

We found similar solutions in the BSBM and BD theories because the dynamics in 
the Einstein frame is very similar in the absence of non-relativistic matter. Likewise, 
one can find identical solutions for the covariant varying speed of light (VSL) theories 
described in Refs. |2niinn!- This does not imply that the VSL, Brans-Dicke and BSBM 
theories are equivalent; merely that one needs to add more general matter (even if 
only as "test matter" in a radiation-dominated universe) for their differences to become 
obvious. Specifically, the coupling to Cem in Q is replaced in Brans-Dicke theory by 
non-minimally coupling the matter fields to the metric. 

In summary, we have considered some simple exactly soluble models for closed 
bouncing universes in theories with varying a and varying G and examined the effects 
of simple non-equilibrium behaviour. Even though the expansion scale factor of the 
universes undergoes periodic oscillations about a finite non-singular expansion minimum, 
we find steady monotonic change in the values of a and G from cycle to cycle both in 
the presence and absence of non-equilibrium behaviour. In the non-equilibrium case the 
oscillations of the expansion scale factor grow monotonically in amplitude and period 
in accord with the Second Law, and the expansion dynamics approaches flatness. 
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